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Oh. We study the sine-Gordon kink diffusion at finite temperature in the overdamped limit. By 

means of a general perturbative approach, we calculate the first- and second-order (in temperature) 
contributions to the diffusion coefficient. We compare our analytical predictions with numerical 
simulations. The good agreement allows us to conclude that, up to temperatures where kink-antikink 
nucleation processes cannot be neglected, a diffusion constant linear and quadratic in temperature 
gives a very accurate description of the diffusive motion of the kink. The quadratic temperature 
r& • dependence is shown to stem from the interaction with the phonons. In addition, we calculate and 

compute the average value {(j>{x,i)) of the wave function as a function of time and show that its 
width grows with ^Jt. We discuss the interpretation of this finding and show that it arises from the 
dispersion of the kink center positions of individual realizations which all keep their width. 
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I. INTRODUCTION 



There is no longer any controversy about the physical relevance of noise effects in spatially extended, nonlinear 
systems Indeed, the pervasive, joint role of nonlinearity and (static or dynamic) disorder has already been 

£j ■ recognized in biophysics, electronics, optics, fluids, condensed matter, computational physics, etc. In most of these 
fields, nonlinear phenomena involve nonlinear coherent excitations, such as solitons or solitary waves, which play a 
key part in the corresponding system dynamics. It is because of this nowadays well established fact that much effort 
has been devoted to understand how stochastic perturbations affect solitons, mostly during the decade of the 80's 
(see for reviews). In fact, early numerical simulations Q already revealed that </> 4 solitary waves underwent 

Brownian-like motion in the presence of additive white noise, i.e., of thermal fluctuations. Subsequent works focused 
£NJ , on the study of soliton diffusion, since it may be crucial in a number of problems, such as photoexcitation dynamics, 
photoconductivity of conducting polymers, or transport by phase solitons in charge-density-wave systems, to name a 
few§. 

Among the different soliton-bearing nonlinear models which have been studied in the above context, one which has 
received a great deal of attention is the sG equation. The interest in this model is both theoretical, as it displays the 
main features of more realistic and complicated cases while remaining analytically tractable, and applied, as it very 

S approximately describes the dynamics of many physically relevant systems, such as one-dimensional (ID) magnets Q 
or long Josephson junctions ||, for instance. Soliton diffusion governed by the sG (and other nonlinear Klein-Gordon 
equations) has been studied along two main, different lines which are discussed and compared, e.g., in Ref. p0[ . 
The first one consists of considering extended excitations of the system (phonons) in equilibrium with both a single 
sG soliton and a heat bath at temperature T . This approach leads to two distinct diffusion regimes: anomalous 
diffusion, characterized by a diffusion constant proportional to T 2 , and viscous diffusion, when the appearance of a 
J> . dynamical damping coefficient yields a diffusion constant proportional to T -1 . We will not follow this approach here; 
the interested reader is referred to the detailed review by Y. Wada pd| ]. The second manner is a la Langevin, i.e., 
introducing the influence of an external thermal bath by means of local fluctuations of the string and a local damping 
force related to that by an appropriate fluctuation-dissipation relationship. The corresponding equation of motion is 
then 

(f>tt ~ <Pxx + sm(<p) = -oujn + r)(x,t), (1) 

with 

(n(x,t))=0 ) (2a) 
(r](x, t)ri(x',t')) = D5(x ~ x')S(t - t'), (2b) 
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where the diffusion coefficient D = 2aki,T, kf, being the Boltzmann constant, and —a<p t being the damping term with 
a dissipation coefficient a. This equation has been considered a number of times in the literature (see, e.g., [|| and 
references therein; see also for related experimental work). 

In this work, we focus on the Langevin version of the problem, with the aim of improving the analytical results 
obtained in the aforementioned works as well as of verifying them by numerical simulations specifically planned to 
that end. Furthermore, we concern ourselves with the overdamped limit of the sG equation, which reads 

at(j)t - (pxx + sin(0) = e?7(x, t, <£,...). (3) 

Note that we have introduced a factor e in front of the noise term for convenience in the analytical calculations in 
section II. This equation (without noise, e = 0) was already considered by Eilenberger in |l3| , as the limit of the 
sG equation ([!]) in the case when the dissipation effect is strong enough in Eq. (|l|) and there is an input of energy 
into the system (see, e.g., JlJ|,[lf| and references therein). On the other hand, Eq. (||), with additive noise as in (0), 
is interesting in itself: For example, it has been proposed as a model for crystal growth (see jl6|-|Tg|l and references 
therein) . Equation (||) has also been studied to analyze the kink contribution to transport properties when the system 
is driven and thermally activated Jl^Jl9| -pT|| . In particular, the work of Kaup J2lj is the most closely related to the 
present one, as it presents a singular perturbation theory to compute the first-order (in T) correction to the kink 
mobility as well as the change of its shape. However, to our knowledge the free diffusion problem for the overdamped 
sG equation has not been adressed in the literature to date and, therefore, we believe that our results will be interesting 
by themselves. On the other hand, we also hope that what we learn in this case can be used towards obtaining a 
more complete, accurate picture of the full sG problem; we will discuss this question in the conclusions. 

The outline of the paper is as follows: In section II, using a general perturbative method |l3| which we recall in 
detail, we calculate the correlation functions of the position and the velocity of the kink center up to second-order in 
kf,T, as well as the diffusion coefficient and the mean value (4>(x, t)) for fixed t. In section III we numerically integrate 
the stochastic partial differential equation (||), with noise given by Eq. (||), using the Heun scheme [^2j and compute 
the time correlation function of the position of the kink center and the diffusion coefficient. We compare these results 
with the theoretical ones obtained in section II and find an excellent agreement. Finally, in section IV we discuss our 
results, summarize our main conclusions, and sketch lines for future research. 



II. A GENERAL PERTURBATIVE APPROACH 

Following the Ansatz proposed in Jl^,|23| , we assume that the solution of Eq. (||) can be expanded as 

/+oo 
dkA k (t)f k [x-X(t)}, (4) 
-oc 

where f k [x — X(t)] are the eigenfunctions of the linearized version of Eq. (0) [with e — 0], which along with Jt[% — 
d(f) 

X(t)] = — — [x — X(t)], form a complete set of orthogonal eigenfunctions (see Appendix I). The first term in the 



expansion (|4|) represents the translational mode related to the position of the kink center X(t), whereas the second 
one characterizes the phonon modes (linear excitations around a kink) of the system. We will focus on the kink center 
motion as described by X(t), as it is well established that such a particle-like picture is very generally enough to 
describe the behavior of the kink as a whole (X playing the role of a collective coordinate; see, e.g., |^4| for a review). 

In order to calculate the dynamics of the kink center, we begin by inserting (||) in (||), and projecting on the 
orthogonal basis [see Appendix I, relationships ([36])] we obtain a system of differential equations for the unknown 
functions X(t) and A k (t): 



X{t) = --X(t) dkA k (t)h{k) - — dk dk'A k (t)A k ,(t)R 3 (k,k r ) + 

' s J-oo 16a J-oo J-oo 

f T [x-X{t)] r)(x,t) dx- 

— QO 

Kj Z'+OC /- + 00 

dk / dk! i dk 2 A k (t)A kl (t)A k2 (t)R 6 (k,k 1 ,k 2 ), (5) 

J —oo J —oo 



/D 
1 

48^ 



8A uj 2 r+oo -i r+oc 

-JL + -±A k (t)=X(t) dkA k (t)I 3 {k,k') + — I dk I dk'A k {t)A k ,{i)R A {k,k') 



2 



■s/D r +oc 

~— fp[x-X(t)} 7](x,t) dx + 

" J-oo 



+00 z' + oo 



where 



/ (9 f 



y— — , /7rfc\ 
'2-7T cosh (^-^"J 



f+oo an ■( 2 2 \2 

/?,</.■./.') : / h{0)^f k {6)r kl {e)d9 = L* h± Mr: A* = *'-*, 

1 AiUkOJk' sinh ( — - 



R 4 (k,k') = [f*,(e)} 2 -£f k (9)d6, Ri(k,k) 



8V27T cosh 



(?) 



R 7 {kMM)= cos&o) r k ,(0)f k (e)r ki (6)f k2 (e)de. (7) 
«/— oc 

We now recall that, if we set e = in (^), the static kink is an exact solution; hence, in what follows we will 
consider e as a small perturbative parameter, and expand A k (t) and X(t) in powers of e. By substituting the series 
Ak{t) = S^Li e ™^fc(*) an( i X(t) — SnLi en X n (t) i n (H) an d @ we find a set of linear equations for the coefficients 
of these series. We only write down here the systems of equations up to order e 3 , leaving out the cumbersome (albeit 
straightforward) equation for A k (t): 

0(e) 

X 1 {t) = e 1 {t) 1 (ei(t)>=0, (e 1 (t)e 1 (t')) = £^6(t-t>), (8a) 

BA 1 j 2 f (t\ D 

*-(*) + i^4(t) = (e fc (t))=0, (e k (t)e k ,(t')) = ^S(t-t')6(k-k'); (8b) 



Of 



0(e 2 ) 



Mt) = dkAl(t)h(k) 



+ OC 



1 

~ 16a" 

ba 2 up- r +0 ° 
—±(t) + -±Al(t)=X 1 (t) / dkAi(t)h(k,k) 
at a /_„ 



rffc / dfc'^(t)A^(t)i? 3 (fc,A;'), (9a) 



-i /- + 00 /-+0O 

— / dfc / *'4(i)4,(i)i? 4 (fc,fc'); (9b) 



0(e 3 ) 



^3(*) = -^# l +C ° dkAl{t)h(k)-^ f + °° dkA\{t)h{k) 



' — OO 

1 r+oo r+oo 

— / dfc / &'A 2 w4,(t)i? 3 (fc,fc') 

16a J^oo J^oo 
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X. 

+00 



dk / dk'Al(t)Al,(t)R 3 (k,k') 
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/ dk I dk\ 
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dk 2 Al{t)Al i {t)A 1 k2 (t)Rs{k,k 1 ,k 2 ). 



We now proceed with the first-order equations. The solutions of fl8a|) and ( |8b|) can be written as 
Xx(t) = /* ei (r)dr, 4(t) = exp(-^) /" exp(^)e fc (r)dr, 
respectively. From these relations we can immediately compute averages over the quantities of interest, such as 



(X 1 (t))=0, {Xi(t)Xi(t')) = — M, 



{x 1 (t))=o, (Xi(t)Xi(t')) = 



8a 2 
D 

8a 2 



(Afc(t)) =0, (Al(t)Al(t>)) 



D 



2auo 2 k 



exp I 



wf(i + t') 



(10) 

(11) 

(12) 
(13) 
(14) 



where M = min(t, t'). For the next orders, the calculations are more involved but not difficult. After some tedious 
algebra, from Eqs. (|9a|)-(|l0|) we find the average values of the position and velocity of the kink center 



(X 2 (t))=0, (X 2 (t))=0, 
(X 3 (t)) = 0, (X 3 (t)) = 0; 



whereas it can be shown that, for large enough times, 

Zak^T 



{\A 2 k {t)\) 



16V2^w 2 



dk 



(nk\ 
' Wfecosh^— j 

The corresponding correlation functions for X 2 (t) and X 2 (t) are 



1.62386. 



(15) 
(16) 



(17a) 
(17b) 



{X 2 (t)X 2 (t')) = 



D 2 M D 2 tt 



512a 3 4096a 2 



exp 



( ~ 2a, 2 M/c 



a - 1 



dk 



u\ cosh 2 (^) 



(X 2 (t)X 2 (t')) = (X^X^t')) x 

Dn r+00 exp ( - uj 2 k \t' - t|/a) - exp ( - u 2 (t' + i)/a) 



256a 



dk 



cosh^ ( — 



(18) 



(19) 



Notice that the cross correlation function of X\(t) and Xs(t) is of the same order as (X^ififX^t')) , and also that 
{X 1 {t)X 2 {t')) = 0. So, from Eqs. (fal) and © we have 



(X 3 (t)X 1 (i'))=(X 2 (t)X 2 (i')) 



I? 2 



256a 3 



dfc Ii(k)x 



{( 



+°° , fl 4 (m,m)\ rM a(exp(-w?M/a) - 1) 
am » — 



(20) 



dn- 



Ri{n,n) a r (exp(-2u 2 M/a) - 1) (exp(-u; 2 M/a) - 1) 



2co 2 
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(X 3 (t)X 1 (t')) = --(X 1 (t)X 1 (t'))x 

f +x r 1 i ( 21 ) 

Finally, from (pi), @ and @ we obtain the final result, namely that for large t [i.e., taking the limit t — ► oo in 
Eqs. ( p^[ ) and ( j2l| ) in all terms except those related to -Xi(t)] the correlation function (X{t)X{t')) is given up to order 
e 4 by 

(X(t)X(t'))=e 2 (X 1 (t)X 1 (t / )) + 

+ eH(X 2 (t)X 2 (t')) + (X^tJXsit')) + (X^X^t'))) + ... 

= + e 2 (^ + Y^ 2 )^} + o(e 4 ). (22) 

We now return to our original equation notation: We set e equal to one and consider y/kbT as the small parameter. 
When t goes to infinity and imposing e = 1, from Eqs. (|l2|), ( p8| ) and (|2^) we find 

Note that the slope of this function is the kink diffusion coefficient, so if one takes into account the second-order 
correction one obtains that the diffusion coefficient is a quadratic function of the temperature. We postpone our 
comments to Sec. IV, where a comparison with the previously available results will be made. 

To complete this work, we can calculate in a very simple way the average value of the wave function <p(x, t) in first 
order: From Eq. (4) we have that 

(<f>(x, t)) = (Mx - eXi(t)]) + 0(e 2 ). (24) 

In this last relation we have taken into account that (Ak(t)) = e{A\(t)) + 0(e 2 ) and (A\(t)) = [see Eq. (p^)]. 

If we solve the corresponding Fokker-Planck equation for X\ [see Eq. (§ah], we obtain that the probability distri- 
bution function for X\ is a Gaussian function given by 




So, one can define the average value (<p(x, t)) as 

/ + OC 
dXip(Xi)c(>o[x - eXxit)]. (26) 
-oo 

Unfortunately we have not found the analytical expression for this integral. But we have calculated it numerically, 
and below we will compare it to the simulations for the full partial differential equation. 



III. NUMERICAL SIMULATIONS 

For our numerical simulations of the partial differential equation (|^), we have used the Heun method [p2| , whose 
stochastic properties are well known and suitable for comparison to our theoretical predictions. We numerically 
integrate Eqs. (|^), with white noise (||), starting from an unperturbed kink at rest and taking the average values over 
1000 realizations. The other parameters are a = 1, Ax = 0.05, and At = 0.001. In the evaluation of the simulations, 
we have defined the center of the kink as follows: We first find all the lattice points % such that 4>i < ir and 4>i+i > tt 
or vice versa. We then interpolate to obtain the points xt where the field cf> crosses tt. In case that, due to the noise- 
induced deformation of the kink, there is more than one such Xi, we average them to finally obtain the numerical 
kink center position, x c . As discussed below, this introduces some error, but other alternatives we tested (such as the 
center of mass, for instance) gave results which did not really represent the kink location, and moreover its calculation 
from numerics is much less accurate. Once the center is obtained, we computed also its variance ([X(t)] 2 ). 
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FIG. 1. Simulations with initial condition given by a static kink initially located at X(0) = 0, and subject to a thermal bath. 
As a continuous (but wiggly) line, we have plotted {[X (t)] 2 ) as obtained by numerical integration of Eq. (^|) for (a) kbT = 0.2, 
(b) ki,T = 0.4, (c) ktT — 0.6 and (d) kbT = 0.8. Overimposed to these lines the linear regression of the numerical results for 
t > 30 is shown (long-dashed line). The solid line is the theoretical prediction ([X(t)] 2 } from Eq. (p3|); this line practically 
overlaps with the linear regression in Fig. a, b, and c. The first-order result {[X(t)] 2 } from Eq. @ is shown as a dot-dash line. 

Figures 0(a)- (d) show a comparison of our numerical results with the analytical predictions, Eqs. ( |l2|) and (53), 
for different values of kbT. We see that there is an excellent agreement between theory and numerics except for the 
highest value of kbT [Fig. 0(d)] . We have checked that this disagreement arises from the way we compute the kink 
center: For such large values of the noise, points where <j)(x, t) = tt are found all over the system, irrespective of their 
distance to the kink center (we note, however, that the temperature was not high enough to create new kink-antikink 
pairs). Those points contribute to the center position through our averaging procedure, and in fact their contribution 
can be shown to be additive, i.e., it amounts to move the whole curve ([X(i)] 2 ) upwards. This is indeed what occurs 
in Fig. |l|(d), and as we will see below the slope is very close to the predicted one. The same behavior is found for 
higher temperatures in so far no new kinks are created (not shown) . Interestingly, a first conclusion that can be drawn 
from these figures is that already for not so high temperatures, kbT > 0.4, as time passes the kink behavior becomes 
more and more different from the first-order prediction, showing clearly the necessity for the second-order correction. 



G 




FIG. 2. Lower solid line: the function D\\ upper solid line: D2, which represent the first- and second-order results for the 
kink diffusion coefficient [see Eqs. ( [L?] ) and (0)]. Diamonds represent the numerical values of the kink diffusion coefficient, 
obtained by numerical integration of Eq. (^) with final time tf = 200 (as in Fig. ^) and different values of kbT. A quadratic 
regression of these numerical values is also plotted (long-dashed line). 



We have calculated the numerical values of the diffusion coefficient for several temperatures by taking the slope of 
([X(t)] 2 ), which we obtain from a linear fit of the data for not so early times (t > 30) to avoid transient effects coming 
from the adjustment of the kink to the heat bath. Note also that our prediction for the second-order contribution 
was obtained in the large-time limit, so we should not try to fit the whole evolution. The figures also show those 
linear regressions. Subsequently, in Fig. we have compared the computed slopes with the first and second-order 

coefficients D\ — — — , and D 2 = — — \ I + ( — + -—v 2 ) kbT > [see. Eqs. (fL3) and (E3)]. The comparison is once 

4a 4a I V32 128 / > 

again very good, and points out very clearly that for values of kbT as low as 0.3, the first-order prediction begins 

to deviate from the diffusion constant measured in the simulations. In addition, the quadratic fit to the simulation 

results, shown as a long-dashed line in Fig. |2[ practically coincides with the second-order prediction in the whole 

studied range. 

As a final verification of our results, in Fig. [| we have plotted the mean value (<f>(x, t)) of the wave function at three 
different times along its evolution, both as obtained from the numerical simulation of the partial differential equation 
and from the numerical evaluation of Eq. (|26|). The perfect agreement between these expressions provides us with 
a hint to derive an approximate analytical estimate of the evolution of (4>) from the integral (pfj|). From Fig. |^ one 
immediately concludes that the width of ((f>) increases with temperature and time. Let us define the width of ((f>) by 
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FIG. 3. Solid lines: Snapshots of the evolution of (<f)(x,t)), obtained from numerical simulations of the partial differential 
equation, for fixed times 40, 120 and 200, respectively. The initial kink (unperturbed, at rest, is also included for comparison) . 
The width of (</>) increases as time progresses. The overimposed points have been computed numerically from the integral (Eq). 
Plots correspond to ktT = 0.4 (a), and 0.8 (b); the width of (0) is seen to increase also with temperature. 



L(t) 



l {[<j> x {x,t)Y)dx 



(27) 



With this definition, we can now calculate {[<f> x {x,t)] 2 ) by using the distribution function of Xi(t); this procedure 
yields 



L{t) 



<[*i(*)] 2 



(28) 



where L\ 



J^dx [* 2 /cosh 2 (x)] 



= 0.8225. 



J^dx [l/cosh 2 (x)] 

It is important to note that, of course, we could define L(t) using ((f> x (x,t)) instead of {[4> x (x, t)] 2 } in the above 
expression, or equivalently another quantity which is localized around the kink center. However, as all possible (and 
sensible) definitions of L(t) give more or less the same results, the difference between them becomes a constant factor 
when ([JTi(i)] 2 ) increases above L\ (for example, for large enough t). So, we expect that the ratio 



L(t) 
L{tf ix ) 



(29) 



for large enough t and tfi x . 

Figure ^ shows a comparison of this prediction with the numerical evaluation of the width of {(f)) from Eq. (|2£ 
From these plots we see that the broadening of (</>) behaves indeed as \ft: We can compare the analytical slope equal 
to 0.5 with the numerical ones equal to 0.4276 and 0.4517 for plots (a) and (b) respectively. The slope in (b) is closer 
to the analytical value due to the fact that tfi X is larger than in case (a), which agrees with the above considerations. 
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FIG. 4. Solid lines: Analytical values of ln[L(t)/L(tfi x )] for tnx = 40 (a) and tfi x = 80 (b). In both cases a = 1 and 
ktT — 0.6. Long-dashed lines: numerical values, calculated from (|26[). The solid lines over the long-dashed ones correspond to 
the linear regression of the numerical points. 



IV. DISCUSSION AND CONCLUSIONS 



As we have seen in the previous section, our second-order theoretical predictions constitute a very accurate descrip- 
tion of the kink dynamics for a wide range of temperatures, up to a value of k^T ~ 1. In fact, the range of validity of 
the analytical results might be somewhat higher, provided a better way to estimate the kink center from the numerical 
simulations could be devised. In any event, the occurrence of 7r-crossings far away from the kink center for values 
around k^T ~ 1 indicates that further increments of the temperature would undoubtedly produce kink-antikink pairs, 
thus invalidating our collective coordinate approach which necessarily relies on the identification of the individual kink 
propagation. We note that this value is a little over 10% of the kink rest mass (Mo = 8 in our units); in this respect, 
a similar result was obtained in |25fl for the overdamped 4 model by means of a similar perturbative approach (with 
the caveat that the numerical data presented in pq| only allow one to guess what is the range of validity of their 
results). 

It is interesting to pursue further the comparison of the results for the sG and cases. In our calculations for 
sG, we have found that the second-order correction is clearly smaller (albeit relevant) than the first-order one. The 
structure of the perturbative calculation allows to identify the origin of that correction: It comes from the interaction 
of the phonons (described by the functions Ak) with the kink. Now, in the </> 4 case, the situation is quite different: 
Indeed, the second-order correction is much larger than the one we find here, and the reason is the so-called internal 
mode, present for </> 4 kinks and absent in the sG case. The coupling between this internal mode (which has been shown 
to act as a reservoir of energy available for exchange with the kink translation mode p6|]) and the kink motion can be 
shown, by a careful examination of the calculation in [§5| , to be responsible for most of the second-order correction, 
while the phonons produce a second-order term comparable to the one we have found. We thus see that, while the 
range of validity of the analytical approach is in principle the same in both cases, the physics is certainly different, 
and in fact the question arises as to the validity of this kind of perturbative calculation for the <j> problem in view 
of the large contribution of the internal mode. This is an interesting question that deserves further analytical and 
numerical work. 

Coming back to our results for the sG kink, the fact that the second-order correction is smaller than the first- 
order term makes us confident that our expansion is likely to be free of problems coming from secular terms. This 
belief is reinforced by the result that, up to the validity range discussed above and limited by kink-antikink creation 
phenomena, the second-order result describes very accurately the kink behavior, which deviates very little from the 
predicted diffusive motion. It is then reasonable to expect higher-order contributions (whose calculation is extremely 
cumbersome, but feasible in principle) to be negligible, thus yielding our theoretical result as the final one for the kink 







diffusion in the overdamped sG problem. In this context, it is also important to realize that Eqs. (pa)-(8b), which are 
only first-order, can also be obtained following the McLaughlin and Scott procedure (see also fl24f). However, the 
advantages of the perturbative scheme we have used are, on one hand, that we were able to obtain the next order in 
the expansion, and on the other hand, we demonstrated that the second order originates in the interaction between 
phonon and translational modes of the sG kink. 

A final remark on our results relates to the mean value of the wave function (4>(x, t)} as a function of t, that must 
not be interpreted as the shape of the kink; in contrast to the interpretation in |2q | . We first note that the width of 
the kink cannot increase from its value when unperturbed at rest; the sG equation, being Lorentz invariant, implies 
that the kink width diminishes when in motion, and therefore an increasing of the width would be very difficult to 
understand on physical grounds. Indeed that is not the case. The broadening of the mean wave function comes in fact 
from the dispersion of the individual realizations, as is immediately seen from Fig. ^| As may be seen all individual 
realizations show a width comparable to the initial kink width, which agrees with our physical intuition. The observed 
\ft behavior, discussed at the end of the preceding section, is then evidently related to the fact that the variance of the 
kink position has that behavior too. The correct interpretation of the width of (cf>(x,t)} is that it represents the area 
in which the kink can be located as its diffusive motion progresses. A similar result has been found for multiplicative 
noise in |28] (see also M and references therein). 




FIG. 5. Average of the wave function for k^T = 0.4 and t = 200 obtained from 1000 realizations (wider solid line), compared 
to the average of only 5 realizations (dot-dashed line) . Also represented are 3 of these individual realizations. Note the different 
slope and width of the average values as compared to individual realizations. 
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To conclude, we want to stress that our main result is the quadratic dependence of the diffusion constant on the 
temperature, stemming from the kink-phonon interactions. This has been verified numerically to a high degree of 
accuracy. We have carried out standard Langevin dynamics simulations following a well grounded procedure, the 
Heun method, as far as statistical properties are concerned p2| . We can thus be sure that what we are dealing with 
is indeed the dynamics of a sG kink at finite temperature. Therefore, our analytical calculations and our numerical 
simulations establish firmly the quadratic dependence of the kink diffusion constant on the temperature for the first 
time. Now the question remains as to the behavior of underdamped sG kinks. Preliminary calculations p^] seem 
to indicate that for underdamped sG kinks the second-order correction is of the same order as that found here, 
which would support the applicability of the previous calculations at least for small temperatures and not too small 
damping. To date, no detailed comparison with numerical simulations has ever been done to check the importance of 
the second-order correction. On the other hand, it would be interesting to compare the results of our approach with 
the theoretical analysis and experiments in |l2| . Such comparison would provide much insight into the importance of 
second and higher-order corrections in actual physical systems. Work along these lines is in progress ]29| | 
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Note added in proof: After acceptance of this paper, we have implemented an improved algorithm for detecting 
the kink center in our code. With this new procedure, no spurious contributions (see discussion below Fig. |l|) to 
the variance appear: Specifically, Fig. [l](d) is largely improved, and the numerical results overlap the theoretical 
prediction, thus confirming the interpretation we have provided of the discrepancy. A detailed report will be given in 

APPENDIX I 

One class of solutions of (||) [with e = 0] is represented by a static kink 

4>o(x, t) = 4 arctan[exp(x)]. (30) 
The perturbations over this equation may be treated by assuming that the solution of (|^) [with e = 0] has the form 

(t>(x,t) =(f)o(x)+ij(x,t), tp(x,t) -C (t>o{x). (31) 
If we substitute Eq. (|3l]) in (^) [with e = 0] and linearize around (j>o{x), we obtain the following equation for ip(x, t) 



dipt = lf>x 



cosh (a;) 



tl>. (32) 



Then, the solution of ( |32| ) may be written as ip(x,t) = (ir) exp ^ — ^, where fk(x) satisfies the eigenvalue 

problem given by 



d 2 f k 



1 



cosh 2 (x) 



dx 2 

This equation admits the following eigenfunctions with their respective eigenvalues 



fk = co 2 k f k . (33) 



M*) = —. T7-v 4 = 0, (34) 
cosn(x) 

expjikx) [fc-Htanh(s)] 2 2 
f k [x) = -= , uj k = 1 + k . (35) 

V 27T UJ k 

Notice, that fr{x) and fk(x) form a complete set of functions with the orthogonality relations 



11 



/+oo p + oo 

fi{x) dx = 8, / f T (x)f k (x) dx = 
-OO J — OO 



r. 

J —oo 



fk(x)f* k ,(x) dx = S(k-k'). 



(36a) 
(36b) 
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